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The single-spin memory effect is considered within a minimal polaron model describing a single- 
level quantum dot interacting with a vibron and weakly coupled to ferromagnetic leads. We show 
that in the case of strong electron- vibron and Coulomb interactions the rate of spontaneous quantum 
switching between two spin states is suppressed at zero bias voltage, but can be tuned through a 
wide range of finite switching timescales upon changing the bias. We further find that such junctions 
exhibit hysteretic behavior enabling controlled switching of a spin state. Spin lifetime, current and 
spin polarization are calculated as a function of the bias voltage by the master equation method. 
We also propose to use a third tunneling contact to control and readout the spin state. 



One of the most promising directions in the fields of 
molecular electronics and spintronics is the experimen- 
tal and theoretical investigation of spin manipulation in 
quantum dots and single molecules. In particular, new 
methods have been recently developed to investigate spin 
states of single atoms and molecules using spin-polarized 
scanning tunneling spectroscopy [l], 0- Motivated by 
such achievements the promising question arises whether 
a single-spin memory effect (including bistability and 
controlled switching between spin states) is possible. 

One of the ways for single spin manipulation is based 
on the interplay of charge, spin, and vibron degrees of 
freedom in molecular junctions. In various experiments 
the signatures of the electron-vibron (c-v) interaction 
have been observed in atomic scale structures 0, 0, S S] • 
In the case of strong e-v interaction the formation of a lo- 
cal polaron can lead to a charge-memory effect, which was 
first predicted long time ago Ml, and has been recently 
reconsidered in more detail [1,0, HI 11, 12, l]|. Neutral 
and charged (polaron) states correspond to different local 
minima of an effective energy surface and are metastablc 
if the e-v interaction is strong enough. By applying 
an external voltage, one can change the charge state of 
this bistable system, an effect that is accompanied by 
hysteretic charge-voltage and current-voltage curves. A 
similar memory effect was found in recent STM experi- 
ments 14 , HI as a multistability of neutral and charged 
states of single metallic atoms coupled to a metallic sub- 
strate through a thin insulating ionic film; the corre- 
sponding polaron model was discussed in 13, Hj]. 

In this Letter we propose an approach to observe 
a single-spin memory effect by combining the polaron 
memory mechanism and the spin-dependent tunneling. 
To this end we consider a single-level and single-vibron 
quantum system between magnetic leads (Fig. [I]). We 
study the case of a symmetric junction with anti-parallel 
magnetizations of left and right leads, besides the third 
electrode can be used as a gate or to probe the spin state. 
The problem can be solved with well controlled approxi- 
mation in the limit of weak coupling to the leads, where 
the master equation for sequential tunneling can be used. 
Thus we focus our major discussion on this limit. 



The Hamiltonian of the single-level polaron (Anderson- 
Holstcin) model is 

H = iud^du + loqo 1 a + A (a^ + a) h + U n^fii 

a 

+ ^2 [i e ika + eipi)c\ k(T Ci k(7 + (Vik<jc\ ka d a + h.c)j . (1) 



Here the first line describes the free electron states with 
energies e CT , the free vibron of frequency u>o, the electron- 
vibron and Coulomb interactions with coupling strength 
A and U, respectively; a is the spin index and h a = d^dg-, 
n — n| + n j . The other terms are the Hamiltonian of the 
leads and the tunneling coupling (i = L, R is the lead 
index, k labels the electronic states) . The bias voltage V 
is introduced through the left and right electrical poten- 
tials, V = (fiL — <Pr- The energy e CT = e a + ei^o includes 
the bare level energies (ej = ej. = e below) and the elec- 
trical potential (fo describing the shift of the central level 
by the gate voltage Vg and the bias voltage drop between 
the left and right lead: <po = ipn + rj((pL — <Pr) + iVgj 
where < ?/ < 1 describes the symmetry of the voltage 
drop across the junction, i] = 0.5 stands for the symmet- 
ric case considered below. 




Figure 1: (Color online) Schematic picture of the considered 
system: a gated single-level quantum dot interacting with a 
vibron and coupled to ferromagnetic leads. 



2 



The coupling to the leads is characterized by the level- 
width function 



T lr7 (e) = 2nY\V tka \ 2 6{e-e lka ). 
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In the wide-band limit considered below, the spin- 
dependent densities of states in the leads and the 
tunneling matrix elements are assumed to be energy- 
independent, so that Tl g and Y R a are constants. The full 
level broadening is given by the sum T a = Tlv + ^Ra- 
Below we consider a symmetric junction with antipar- 
allel magnetization of the leads and use the notation 
Tli = Tri = r for majority spins and = Tr^ = kT 
for minority spins, k«1. 

The spin effects addressed are particularly pronounced 
in the limit U — > oo, i.e. we neglect the doubly occu- 
pied state, so that only three states in the charge sector 
should be considered: neutral |0), charged spin-up | f) 
and charged spin-down | J.). Using the polaron (Lang- 
Firsov) [3, E3, EH canonical transformation, the eigen- 
states of the isolated system (T = 0) are 
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where the quantum number q characterizes vibronic 
cigenstatcs, which are superpositions of states with dif- 
ferent number of bare vibrons. 

Taking into account all possible single-electron tunnel- 
ing processes for both leads, we obtain the incoming and 
outgoing tunneling rates 



r ?9' ~~ X r "79' ~~ X] Ticr \ M iq'\ 2 fi{ E aq - E 0q') 
i=L,R i=L,R 

= T^lM^ff^^+Mq-q')), (6) 

i=L,R 

r qq' = X = X Tia \ M ii'\ 2 i 1 ~ fi( E °q' ~ E 0q)) 

i=L,R i=L,R 

= J2 F - I^VI 2 (1 - f?& "o(q q'))) ■ (7) 

i=L,R 

Here (e) is the equilibrium Fermi function in the lead 
shifted by the external potential, ff(e) = f°(e—e<pi), and 
M qq i is the Franck-Condon matrix element that can be 
calculated analytically (see Refs. [H, 0, EH, 0, H3| for 
details of the master equation method and calculation of 
the tunneling rates). The incoming rate T^, describes 
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Figure 2: Inverse lifetime 7 <T °/r of the neutral state (thin 
solid line) and the inverse spin lifetime (r CT r) _1 (thick gray 
solid line) as a function of the scaled electron-vibron coupling 
\/u>o at eo = \ 2 /2ujo and at eo = 0.1A 2 /ojo (corresponding 
dashed lines), T = O.lwo- 



tunneling of one electron with spin a from the lead to 
the dot changing the state of the dot from |0g') to \<rq). 
The outgoing rate V®^, corresponds to the transition from 
\aq') to |0g). 

In the sequential tunneling regime the master equation 
for the probability Pg(t), n = 0, |, J. to hnd the system 
in one of the pola ron eigenstates (3j , Q can be written 
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Here the first term describes the tunneling transition into 
the state \nq) and the second term the transition out of 
the state \nq). I V [P] is the vibron scattering integral 
describing the relaxation of the vibrons to the thermal 
equilibrium. 

Finally, the average charge and the spin polarization 
are 
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respectively, and the average current (from the left or 
right lead) reads 
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To proceed further, we calculate the characteristic life- 
times of the neutral, spin-up, and spin-down ground 
states (q = 0). We define the switching rates 7 ' from 
the neutral to the charged state with spin a and vice- 
versa as the sum of the rates of all possible processes 
which change these states 
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Figure 3: (Color online) Inverse spin lifetime as a function 
of normalized bias voltage eV/uio at \/ujq = 3, k = 0.01, 
eo = A 2 /2tt)o for the spin-up state (thin red solid line) and the 
spin-down state (thick blue solid line) and the same for a less 
polarized junction (k = 0.1, dashed lines). 
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Figure 4: (Color online) Populations of the spin-up state 
(red), spin-down state (blue), and spin polarization (green) 
as a function of normalized voltage eV/uio at X/uuo = 3 
and eo = \ 2 /2ujo, the solid (dashed) for increasing (decreas- 
ing)voltage. Inset: sketch of voltage time-dependence. 



In the sequential tunneling approximation the spin life- 
time t„ is determined by the lifetime of the charged state, 
from (JT]), (fTTj) . It reads (assuming that the Fermi energy 
in the leads is zero, see details in 12j, g = (A/wq) 2 ) 



r- 1 = I aa = (1 + K)T ^ff + ^1) ■ (12) 

9 

At large g the sequential tunneling rates are exponen- 
tially suppressed and the cotunncling contribution to t~ x 
becomes dominant. It can be estimated as 23] 



_ 1(ct) _ kV 2 TuI 
~ A 4 ' 



(13) 



Although the cotunncling contribution is not suppressed 
exponentially by Franck-Condon blockade, it is of the 
second order in the tunneling coupling and suppressed 
additionally by the small polarization parameter k and 
large A. At typical parameters, considered in this Let- 
ter, the cotunneling contribution can be neglected, but it 
can be essential at larger tunneling couplings and larger 
temperatures. 

The dependence of t" 1 and 7 <T ° on the scaled electron- 
vibron interaction constant ^fg = X/ujo is shown in Fig. [21 
For large values of A the tunneling from the neutral state 
to the charged state and vice versa is suppressed com- 
pared to the bare tunneling rate T. Hence all states 
arc (meta)stable at low temperatures and zero voltage. 
Moreover, the lifetime of the charged states can be much 
larger than that of the neutral state. 

Next we address the other important question, whether 
fast switching between the two spin states is possible. 
To this end we consider what happens, if one sweeps the 
voltage with different velocities, r exp is the characteristic 
time of the voltage change. At this point an assump- 
tion about the relaxation time ry of the vibrons without 



change of the charge state is due. We assume that the 
relaxation is fast, Ty -c T ai T exp , so that after an elec- 
tron tunneling event the system relaxes rapidly into the 
vibronic ground state |ct0) or |00). In this case the prob- 
abilities P a = J2 q Pq of the charged and P° = J2 q Pq of 
the neutral state are determined from the equations 



^! = ^( 7 o^_ r o p o ) 



dt 
dP 



dt 



_ 7 0<T P< T ; 



(14) 
(15) 



where the switching rates 7 ' , 7 0<T at finite voltage are 
calculated from Eqs. (|6l7llip : 

q 

+ T Ra f°(i / a +cu q + veV)], (16) 

7 0a = E [riCT/0 + woq + (1 ~ v)eV) 

q 

+ rn,f(-?a + woq-VeVj\. (17) 

The voltage dependence of the inverse spin lifetime is 
depicted in Fig. [31 If the voltage is large enough, the 
Franck-Colomb blockade is overcome and the system is 
switched into spin-up (spin-down) state at positive (nega- 
tive) voltage. If the bias voltage is swept fast enough, i.e. 
faster than the spin lifetime at zero voltage, r exp <JC t ct (0), 
both spin states can be considered as stable at zero volt- 
age and hysteresis takes place. This is shown in Fig. [4] 
where the solid (dashed) lines mark the spin population 
for increasing (decreasing) bias voltage. In the oppo- 
site (adiabatic) limit the voltage change is so slow that 
the system relaxes into the equilibrium state, and the 
population- voltage curve is single- valued (Fig. [5]) . 
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Figure 5: (Color online) Spin polarization as a function of nor- 
malized voltage eV/coo at X/u>o = 3 and eo = A 2 /2o)o for three 
different sweep velocities relative to that in Fig(4](here shown 
by green): faster (black), slower (blue) and in the adiabatic 
limit (red dashed line). 



Finally, we study the signatures of the spin polariza- 
tion in the charge current which is most easily accessible 
to experiments. In Fig. [5] we show the bias current and 
the test current to the additional ferromagnetic electrode, 
very weakly coupled to the system, so that it does not 
perturb the state. At large negative voltage applied to 
the electrode the current is sensitive to the orientation 
of the magnetization in the test electrode, thus the spin 
state can be controlled during the experiment. Also such 
a small current can be used to readout the memory ele- 
ment. 

In conclusion, we considered a single-spin memory ef- 
fect and switching phenomena in the framework of a 
single-level quantum dot polaron model, taking into ac- 
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Figure 6: (Color online) Bias current (top black) and test cur- 
rent for spin-up (middle red) and spin-down (bottom blue) 
test electrode magnetization as a function of normalized volt- 
age eV/ido; all other parameters are the same as in Fig. 2] 
The curves are shifted along vertical axes and the magnitude 
of the test current is much smaller than of the bias current. 



count non-stationary effects, in particular the interplay 
between the timescales of voltage sweeping and the quan- 
tum switching rates of meta-stable states. We showed 
that the bistability arises when the quantum switching 
between two spin states is suppressed due to the Franck- 
Condon blockade. Controlled switching of the spin can 
be achieved by applying finite voltage pulses. 
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Collaborative Research Center SFB 689. 
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